/ BBS in a Box 10 / BBS In A Box Volume X (AMUG) (January 1994).bin / Files / Prog / U-Z / VideoToolBox.cpt / Utilities / Quick3 / LogLikelihood.c

< prev

next >

MacOS 8.1

Win98

DOS

browse contents

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 00 0f 4c 6f 67 4c 69 6b | 65 6c 69 68 6f 6f 64 2e |..LogLik|elihood.|
|00000010| 63 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |c.......|........|
|00000020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 54 45 58 54 4b 41 48 | 4c 01 00 00 00 00 00 00 |.TEXTKAH|L.......|
|00000050| 00 00 00 00 00 03 ce 00 | 00 01 7e a7 af 62 f0 a7 |........|..~..b..|
|00000060| af 62 f0 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.b......|........|
|00000070| 00 00 00 00 00 00 00 00 | 00 00 81 81 c7 13 00 00 |........|........|
|00000080| 2f 2a 0d 4c 6f 67 4c 69 | 6b 65 6c 69 68 6f 6f 64 |/*.LogLi|kelihood|
|00000090| 2e 63 0d 43 6f 70 79 72 | 69 67 68 74 20 28 63 29 |.c.Copyr|ight (c)|
|000000a0| 20 31 39 39 30 20 44 65 | 6e 69 73 20 47 2e 20 50 | 1990 De|nis G. P|
|000000b0| 65 6c 6c 69 0d 34 2f 37 | 2f 39 30 20 0d 43 61 6c |elli.4/7|/90 .Cal|
|000000c0| 63 75 6c 61 74 65 20 6c | 6f 67 20 6c 69 6b 65 6c |culate l|og likel|
|000000d0| 69 68 6f 6f 64 20 6f 66 | 20 74 68 65 20 64 61 74 |ihood of| the dat|
|000000e0| 61 20 66 6f 72 20 61 20 | 67 69 76 65 6e 20 70 73 |a for a |given ps|
|000000f0| 79 63 68 6f 6d 65 74 72 | 69 63 20 66 75 6e 63 74 |ychometr|ic funct|
|00000100| 69 6f 6e 2e 0d 49 66 20 | 74 68 65 20 70 73 79 63 |ion..If |the psyc|
|00000110| 68 6f 6d 65 74 72 69 63 | 20 66 75 6e 63 74 69 6f |hometric| functio|
|00000120| 6e 20 69 6e 64 69 63 61 | 74 65 73 20 74 68 61 74 |n indica|tes that|
|00000130| 20 74 68 65 20 70 61 72 | 61 6d 65 74 65 72 73 20 | the par|ameters |
|00000140| 61 72 65 20 6f 75 74 20 | 6f 66 20 62 6f 75 6e 64 |are out |of bound|
|00000150| 73 2c 0d 74 68 65 6e 20 | 4c 6f 67 4c 69 6b 65 6c |s,.then |LogLikel|
|00000160| 69 68 6f 6f 64 28 29 20 | 72 65 74 75 72 6e 73 20 |ihood() |returns |
|00000170| 74 68 65 20 6c 6f 77 65 | 73 74 20 70 6f 73 73 69 |the lowe|st possi|
|00000180| 62 6c 65 20 6c 6f 67 20 | 6c 69 6b 65 6c 69 68 6f |ble log |likeliho|
|00000190| 6f 64 2c 20 6d 69 6e 75 | 73 20 69 6e 66 69 6e 69 |od, minu|s infini|
|000001a0| 74 79 2e 20 0d 49 66 20 | 79 6f 75 72 20 6d 61 63 |ty. .If |your mac|
|000001b0| 68 69 6e 65 20 64 6f 65 | 73 6e 27 74 20 64 65 61 |hine doe|sn't dea|
|000001c0| 6c 20 69 6e 20 69 6e 66 | 69 6e 69 74 69 65 73 2c |l in inf|inities,|
|000001d0| 20 74 68 65 6e 20 72 65 | 70 6c 61 63 65 20 31 2e | then re|place 1.|
|000001e0| 30 2f 30 2e 30 20 62 79 | 20 44 42 4c 5f 4d 41 58 |0/0.0 by| DBL_MAX|
|000001f0| 20 69 6e 0d 74 68 65 20 | 64 65 66 69 6e 69 74 69 | in.the |definiti|
|00000200| 6f 6e 20 6f 66 20 49 4e | 46 2e 0d 0d 54 68 65 20 |on of IN|F...The |
|00000210| 74 68 69 72 64 20 61 72 | 67 75 6d 65 6e 74 20 69 |third ar|gument i|
|00000220| 73 20 61 20 70 6f 69 6e | 74 65 72 20 74 6f 20 61 |s a poin|ter to a|
|00000230| 20 70 73 79 63 68 6f 6d | 65 74 72 69 63 20 66 75 | psychom|etric fu|
|00000240| 6e 63 74 69 6f 6e 2c 20 | 65 2e 67 2e 20 57 65 69 |nction, |e.g. Wei|
|00000250| 62 75 6c 6c 2e 63 0d 0d | 48 49 53 54 4f 52 59 3a |bull.c..|HISTORY:|
|00000260| 0d 38 2f 32 34 2f 39 31 | 20 64 67 70 20 63 6f 73 |.8/24/91| dgp cos|
|00000270| 6d 65 74 69 63 0d 2a 2f | 0d 23 69 6e 63 6c 75 64 |metic.*/|.#includ|
|00000280| 65 20 22 51 75 69 63 6b | 33 2e 68 22 0d 0d 23 64 |e "Quick|3.h"..#d|
|00000290| 65 66 69 6e 65 20 49 4e | 46 20 28 31 2e 30 2f 30 |efine IN|F (1.0/0|
|000002a0| 2e 30 29 0d 0d 64 6f 75 | 62 6c 65 20 4c 6f 67 4c |.0)..dou|ble LogL|
|000002b0| 69 6b 65 6c 69 68 6f 6f | 64 28 64 61 74 61 52 65 |ikelihoo|d(dataRe|
|000002c0| 63 6f 72 64 20 2a 64 61 | 74 61 2c 70 61 72 61 6d |cord *da|ta,param|
|000002d0| 52 65 63 6f 72 64 20 2a | 70 61 72 61 6d 73 2c 0d |Record *|params,.|
|000002e0| 09 50 73 79 63 68 6f 6d | 65 74 72 69 63 46 75 6e |.Psychom|etricFun|
|000002f0| 63 74 69 6f 6e 50 74 72 | 20 50 73 79 63 68 46 75 |ctionPtr| PsychFu|
|00000300| 6e 29 0d 7b 0d 09 64 6f | 75 62 6c 65 20 4c 4c 2c |n).{..do|uble LL,|
|00000310| 70 3b 0d 09 6c 6f 6e 67 | 20 6e 3b 0d 09 69 6e 74 |p;..long| n;..int|
|00000320| 20 69 3b 0d 09 63 6f 6e | 74 72 61 73 74 52 65 63 | i;..con|trastRec|
|00000330| 6f 72 64 20 2a 63 50 74 | 72 3b 0d 0d 09 4c 4c 3d |ord *cPt|r;...LL=|
|00000340| 30 2e 30 3b 0d 09 66 6f | 72 28 69 3d 30 3b 69 3c |0.0;..fo|r(i=0;i<|
|00000350| 64 61 74 61 2d 3e 63 6f | 6e 74 72 61 73 74 73 3b |data->co|ntrasts;|
|00000360| 69 2b 2b 29 7b 0d 09 09 | 63 50 74 72 3d 26 64 61 |i++){...|cPtr=&da|
|00000370| 74 61 2d 3e 63 5b 69 5d | 3b 0d 09 09 70 20 3d 20 |ta->c[i]|;...p = |
|00000380| 28 2a 50 73 79 63 68 46 | 75 6e 29 28 63 50 74 72 |(*PsychF|un)(cPtr|
|00000390| 2d 3e 63 6f 6e 74 72 61 | 73 74 2c 20 70 61 72 61 |->contra|st, para|
|000003a0| 6d 73 29 3b 0d 09 09 69 | 66 28 70 3d 3d 49 4c 4c |ms);...i|f(p==ILL|
|000003b0| 45 47 41 4c 5f 50 41 52 | 41 4d 45 54 45 52 53 29 |EGAL_PAR|AMETERS)|
|000003c0| 72 65 74 75 72 6e 20 2d | 49 4e 46 3b 0d 09 09 6e |return -|INF;...n|
|000003d0| 3d 63 50 74 72 2d 3e 63 | 6f 72 72 65 63 74 3b 0d |=cPtr->c|orrect;.|
|000003e0| 09 09 69 66 28 6e 3e 30 | 4c 29 20 4c 4c 20 2b 3d |..if(n>0|L) LL +=|
|000003f0| 20 6e 2a 6c 6f 67 28 70 | 29 3b 0d 09 09 6e 3d 63 | n*log(p|);...n=c|
|00000400| 50 74 72 2d 3e 74 72 69 | 61 6c 73 20 2d 20 63 50 |Ptr->tri|als - cP|
|00000410| 74 72 2d 3e 63 6f 72 72 | 65 63 74 3b 0d 09 09 69 |tr->corr|ect;...i|
|00000420| 66 28 6e 3e 30 4c 29 20 | 4c 4c 20 2b 3d 20 6e 2a |f(n>0L) |LL += n*|
|00000430| 6c 6f 67 28 31 2e 30 2d | 70 29 3b 0d 09 7d 0d 09 |log(1.0-|p);..}..|
|00000440| 72 65 74 75 72 6e 20 4c | 4c 3b 0d 7d 0d 0d 00 00 |return L|L;.}....|
|00000450| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000460| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000470| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000480| 00 00 01 00 00 00 01 4c | 00 00 00 4c 00 00 00 32 |.......L|...L...2|
|00000490| 0a 8e 00 00 66 2a 00 00 | 27 fc 00 00 00 00 00 00 |....f*..|'.......|
|000004a0| 00 00 00 08 00 00 00 00 | 00 00 00 0e 00 ac 00 18 |........|........|
|000004b0| 0f 4c 6f 67 4c 69 6b 65 | 6c 69 68 6f 6f 64 2e 63 |.LogLike|lihood.c|
|000004c0| 02 00 00 00 54 45 58 54 | 4b 41 48 4c 00 00 00 00 |....TEXT|KAHL....|
|000004d0| 00 00 54 45 58 54 4b 41 | 48 4c 00 00 00 00 00 00 |..TEXTKA|HL......|
|000004e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000004f0| 00 00 a7 e5 d5 9f 00 00 | 00 00 00 00 01 7e 00 00 |........|.....~..|
|00000500| 01 10 00 00 0b b8 00 00 | 01 74 00 00 00 00 00 00 |........|.t......|
|00000510| 01 50 00 00 06 ec 00 00 | 02 e8 00 00 00 00 00 00 |.P......|........|
|00000520| 00 00 00 28 00 00 00 00 | 00 00 00 46 01 cc 00 3c |...(....|...F...<|
|00000530| 00 00 00 00 01 66 00 03 | 00 09 80 00 00 78 00 00 |.....f..|.....x..|
|00000540| 0a c8 00 00 15 15 00 01 | 00 00 a2 42 da 21 a2 42 |........|...B.!.B|
|00000550| da 21 a2 42 da 21 00 16 | 00 02 00 00 00 62 42 52 |.!.B.!..|.....bBR|
|00000560| 45 4e 00 07 42 52 45 4e | 54 2e 43 00 00 01 a2 42 |EN..BREN|T.C....B|
|00000570| d1 1a a2 42 d1 1a a2 42 | d1 1a 00 14 00 02 00 00 |...B...B|........|
|00000580| 00 00 00 48 00 09 4d 6f | 6e 61 63 6f 00 68 6f 6f |...H..Mo|naco.hoo|
|00000590| 64 2e 63 00 00 00 00 18 | 00 01 a7 a8 00 01 00 00 |d.c.....|........|
|000005a0| 03 e9 00 00 00 00 00 06 | 00 04 00 3c 00 24 01 dd |........|...<.$..|
|000005b0| 02 3d 00 3c 00 24 01 dd | 02 3d a3 52 2b 4a 00 00 |.=.<.$..|.=.R+J..|
|000005c0| 01 25 00 00 01 98 00 00 | 00 00 01 00 00 00 01 00 |.%......|........|
|000005d0| 00 00 01 4c 00 00 00 4c | 00 00 00 32 00 5d 6f 30 |...L...L|...2.]o0|
|000005e0| 09 ec 00 00 00 1c 00 32 | 00 00 4d 50 53 52 00 00 |.......2|..MPSR..|
|000005f0| 00 0a 03 ed ff ff 00 00 | 00 00 00 5f 15 f4 00 00 |........|..._....|
+--------+-------------------------+-------------------------+--------+--------+